o 

(N 



a 
o 

(N 
> 

\o 
q 

o 



X 



On the Efimov Effect in Higher Partial Waves 

K. Helfrich and H.-W. Hammer 

Helmholtz-Institut fiir Strahlen- und Kernphysik (Theorie) 

and Bethe Center for Theoretical Physics, 

Universitdt Bonn, 53115 Bonn, Germany 

(Dated: September 23, 2011) 



Q^' Abstract 

^^ . Using the framework of effective field theory, we present a detailed study of the Efimov effect in 

higher partial waves for systems of two identical particles and a third distinguishable particle. De- 

^\i \ pending on the total angular momentum L, the two identical particles must be bosons or fermions. 

We derive analytical expressions for the elastic and inelastic atom-dimer scattering cross sections 

t^ ■ as well as the atom-dimer relaxation rate at the dimer breakup threshold. For the experimentally 

\^- most relevant case of P-waves, we numerically calculate the atom-dimer scattering cross sections 

and relaxation rates as a function of the scattering length, three-body parameter, and mass ratio 

S \ for energies below breakup threshold. 

a^, PACS numbers: 34.50.-s, 67.85.Pq 



I. INTRODUCTION 

The last few years saw an enormous progress in the field of Efimov physics, theoretically 
as well as experimentally [1[. Already in 1970, Efimov predicted the existence of universal 
three-body bound states with a geometric spectrum for identical bosons at infinite scattering 
length a [2(] 

4")=(,-2./.o)— ;,2^2/^^ (1) 

where m is the mass of the particles, Sq ~ 1.00624 is a transcendental number, and k* is the 
binding wave number of the state labelled n*. This leads to an accumulation point at ii^ = 
and a discrete scale invariance with a scaling factor of e'^'^° ~ 22.7. Away from unitarity, 
there is only a finite number of bound states. As the inverse scattering length is varied 
from negative to positive values, an Efimov state appears at the three-particle scattering 
threshold at a = a_, crosses the unitary limit at 1/a = with the energy given in Eq. ([1]), 
and vanishes at a = a* through the atom-dimer threshold. Whenever an Efimov trimer hits 
one of these thresholds, recombination rates get enhanced |3|-l5|. 

Measuring this enhancement, the Efimov effect has by now been observed in bosonic 
isotopes of cesium |6|,l7|, potassium |8|, and lithium |9|,|l0|. The case of fermionic lithium with 
three different hyperfine states is similar to the bosonic case in many aspects. Because of the 
three different pair scattering lengths involved, however, the phenomenology is much richer. 



After the indirect observation via resonant enhancement of the recombination rates [ll|, |12 
the direct association of Efimov trimers was also achieved |l3l . Il4l . 



Mixed systems also allow for the appearance of the Efimov effect as long as at least 
two of the scattering lengths are resonant. This is most easily realised with two different 
atomic species, one of which has to be bosonic for the S-wave case. If the bosonic species 
is much heavier than the other one, the scaling factor can become considerably smaller and 



thus more favorable for the experimental investigation 15|]. This case has been considered 



in detail in [16[. The first experimental measurement of the Efimov effect in a bosonic 



rubidium-potassium mixture was reported in [17 



However, if the heavy species in a heteronuclear mixture is fermionic, the Efimov effect 
is only present in a P-wave channel and only if the masses differ by at least a factor of 
13.61 [l5|. This behaviour can be generalized to higher partial waves. The parameters s^ 
determining the scaling factor for total angular momentum L can be estimated by [l5| 



sl^sl-LiL + 1), (2) 

where Sq characterizes the corresponding scaling factor for zero angular momentum. s| > 
must be fulfilled for the occurrence of the Efimov effect (also see more detailed discussion 
in section [IT])- The necessary mass ratios for the P-wave case {L = 1) could be realised for 
example with ^Li-*^Sr, '^/^Li-^^^Ba, '^/^Li-^^^Er, or ^/'^Li-^'^^/^'^^Yb. Higher partial waves would 
only become accessible if hydrogen or helium atoms can be used. Dimer-dimer scattering in 
heteronuclear mixtures showing a P-wave Efimov effect was investigated in |l8l ]. 

The main focus of this paper is also on such mixtures, prepared as atom-dimer systems. 
We study in detail how the Efimov effect in higher partial waves affects observables such as 
atom-dimer scattering and atom-dimer relaxation. We briefly introduce our effective fleld 
theory framework and derive analytical expressions for the elastic and inelastic atom-dimer 
scattering cross sections at the dimer breakup threshold. For the P-wave case, we numerically 
calculate the atom-dimer scattering cross sections and relaxation rates as a function of the 



scattering length, three-body parameter, and mass ratio below dimer breakup threshold. 
Finally, some quantities for systems without Efimov effect are also computed. 



II. FRAMEWORK 

In the following, we investigate various heteronuclear atomic systems in detail and closely 



follow the formalism and conventions of Ref. 16 . For convenience, we set ^ = 1 but restore 



the dimensions for final results. We consider systems consisting of two different atomic 
species, where the occurring trimers, Efimov or non-Efimov, are built of one atom of type 
1 and two atoms of type 2. The unlike particles have a resonant S-wave interaction which 
can be tuned using a Feshbach resonance whereas the interaction between identical particles 
can be neglected. The corresponding Lagrangian is given by [16| 



'^■+S*-+*K'''^^'*= 



c = *1 (i 

+ 92dU - 92 {d^^i^2 + *l*|ci) - ^dU^l^2 + ■ ■ ■ , (3) 

where mi/2 denotes the mass of particles of species 1/2, (72/3 are the bare two-body and three- 
body coupling constants, and d is an auxiliary field for a dimer consisting of particle species 
1 and 2. The first two terms in Eq. ([3]) are the kinetic terms for species 1 and 2, while the 
third term is the kinetic term for the non-dynamical auxiliary dimer field which mediates 
the interactions between the particles of species 1 and 2. It becomes dynamical through 
the coupling to particle loops. The fourth and fifth terms generate the S-wave contact 
interactions between particles of species 1 and 2 and a three-body interaction between one 
particle of species 1 and two particles of species 2|j The ellipses stand for higher-order terms 
containing more fields and/or derivatives. As in Ref. |l6|, we denote the ratio of the masses 
of particle species 1 and 2 by 5 = mi/ 1712. In this study, however, we explicitly focus on 
higher partial waves with total angular momentum L > 0. 

Because of symmetry reasons, the Efimov effect can only occur in even angular momentum 
channels if the two like particles are bosons and in odd angular momentum channels if 
they are fermions |l9|. In the following, we will refer to the first case as bosonic and the 
second case as fermionic for simplicity. The nature of the third particle is not relevant for 
our purpose. For inverse mass ratios 6~^ larger than a critical ratio 6~l^, Efimov physics 
can be observed 15|. At S~l^ and beyond, the angular momentum barrier is overcome 
by the attractive interaction between unlike particles (cf. Eq. ([2])). The light particle can 
be thought of as an exchange particle between the two heavy atoms. In this case, the 
"fall to the center" phenomenon typical for the Efimov effect can happen. A (hybrid) 
Born- Oppenheimer description has been used in the limit of a very light particle of species 
1 ISl, l20|. In the case of L = 1, the mass ratio must satisfy 6"^ > 13.61 = 6~l [l5|, |2l|. The 
D-wave Efimov effect starts at 5~^ > 38.63 = 6~2 (l5|, |22| and its observation would always 



be obscured by the already present S-wave effect. Consequently, the fermionic P-wave case 
is the only relevant one in cold atom experiments besides the S-wave bosonic case. The mass 
ratios 6~^ for some possible mixtures showing the P-wave Efimov effect are given in Table [B 



^ To this order, the eocfficient of the dimer kinetic term is not independent of the interaction between 
particles 1 and 2 and can be chosen as 172 for convenience. 



S^Sr i37Ba i67Er i^iyb i^Syb 

^Li 14.5 22.8 27.8 28.5 28.8 

^Li - 19.6 23.9 24.4 24.7 

TABLE I. Mass ratios 6~^ for possible experimental mixtures showing the P-wave Efimov effect. 
It is straightforward to derive Feynman rules and to obtain the full dimer propagator from 



the Lagrangian in Eq. (jS]) (for more details, see Refs. |l6l . l23j ) . In the three-body system, 
we obtain an integral equation for the off-shell atom-dimer scattering amplitude Al{p, k; E) 
known as the STM equation [2^]. The amplitude depends on the relative momenta of the 
atom and the dimer in the initial state, /c, and in the final state, p, as well as on the total 
energy E. All three-body observables can be obtained from this amplitude evaluated in 
appropriate kinematics. Projecting on total angular momentum L, the equation becomes 

AL{p,k;E) = {±l) -i — (-l)^g '^ 



ajj,'^ pk \ 2pkjj,/mi 

TT^ Jq p V ^pqf^/rrii ) 

^ Al{(iME) 

^ -l/a+^-2iJi{E-q'/{2iiAD))' 



(4) 



where the prefactor +1 corresponds to the atoms of species 2 being bosons and —1 to 
fermions, respectively. Moreover, /x = mim2/{mi + 7712) is the reduced mass of two unlike 
atoms, fiAD = m2{mi + m2)/(mi + 2m2) is the reduced mass of an atom and a dimer, 

Ql{z) = 1 Tdx^^ (5) 

2 J^i z — X 

is a Legendre function of the second kind, and Pl{x) is a Legendre polynomial. The log- 
periodic dependence of the three-body interaction g^ on the cutoff has been used to absorb 
the complex three-body parameter into the momentum cutoff A^ in Eq. (jl]) |5|. The inte- 
gration in Eq. (^ is along a straight line in the complex plane from zero to A^. The physical 
values of the amplitude A are obtained by evaluating Eq. (jlj) once more for real p = k 
after the solution along the complex contour has been obtained. The absolute value of the 
cutoff Ac is proportional to the binding momentum of the deepest Efimov state, whereas 
the complex phase determines the width of the Efimov trimers, Ac oc e^^'^^^n^. Physically, 
the parameter i]^ takes into account the effects of deeply-bound dimers which provide decay 
channels for the Efimov trimers |5| . For S-wave Feshbach resonances, r/^, is typically of the 
order 0.1. In the limit of small ?7*/sl, the width of Efimov trimers is given by 



sl V * 2//a^ 



^(5. + ^). (6) 



where Bt is the trimer energy (cf. Eq. (259) in Ref. [23|). Thus for small rf^/si^ the 
width becomes small. Without knowing 77^,, a quantitative prediction is not possible, but we 
generally expect lifetimes comparable to the S-wave case if sl is of order one. 

In the channels without the Efimov effect (odd angular momenta for bosons, even angular 



momenta for fermions, or for 5 ^ < 5 ]), the absolute value of the momentum cutoff Ac can 



be taken to infinity. The three-body interaction in Eq. (|3]) is higher order in these channels 
and all observables are in leading order determined by the scattering length alone. 



III. SCALING FACTOR AND RESONANCE POSITIONS 



If more than one Efimov resonance feature can be measured in an experiment, the scaling 
factor exp^Ti/sL) can be deduced. The quantity sl can be computed analytically by con- 
sidering the large momentum behaviour of Eq. (|1]) |25l-l27l| . In this limit the energies and 
inverse scattering lengths can be neglected compared to the momenta p and q, the inhomo- 
geneous term as well as purely polynomial terms in the integral kernel are suppressed, and 
the momentum integration can be extended to infinity. This leads to the equation 



TT/iV /i io Q \2pqfi/miJ 



° \2pqfi/mi) 



(7) 



where we have defined Al{p) = P-A.l{p ^ k\ E) . Since the equation is scale invariant, it has 
power law solutions. If the Efimov effect is present, the exponent is complex: Al{p) oc 
p±isL _ Identifying the right hand side of Eq. ([7]) as a Mellin transform, we then obtain a 
transcendental equation for sl- 
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where we defined 
and 



1 2 sin [{isL + 2m - L + 2k)(j)] 

ml {L-2k- m)\ isL + 2m-L + 2k cos [{isi + 2m - L + 2A;)f ] 



(p = arcsm 
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if L is even 
if L is odd 



In the case L = 1, Eq. ([8]) reduces to 



sin[(zsi — 1) 



sin[(zsi + 1^ 



2 sin cos [isi — 1 cos[(zsi — l)7r/2] isi + lcos[{isi + l)7r/2] 



(8) 



(9) 



(10) 



(11) 



The corresponding equation for L 



1 



ism 
2- 



(pcos(p 
4/3 sin' 



2 is given by 
sin[(zs2 - 2)(, 



sin[(is2 + 2) 



is2 — 2 cos[(zs2 
sin[2S20] 



2)7r/2] zs2 + 2cos[(is2 + 2)7r/2] 



is2 cos[2S2vr/2] 



(12) 



An equivalent equation for general L using hypergeometric functions was derived by 
Nielsen and coworkers [l9|. The results of Eqs. fITT]) and f lT2|) coincide with the ones obtained 
by making use of Eq. (117) in 19|. The critical mass ratios 6c,l are obtained by observing 

38.63 in agreement with previous 



when Sl tends to zero. We find 5^ I 



13.61 and 5^2 
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FIG. 1. Scaling factors exp(7r/sL) for bosons (b) and fermions (f) as a function of the mass ratio 
6 = mi/m2 for L = 0, 1,2, respectively. The vertical dotted lines indicate the critical mass ratios 
(5c,i and 5c, 2 for the Efimov effect with L = 1,2. 



determinations |2ll . |22| . The resulting scaling factors exp(7r/si) are shown in Fig. [T] as a 
function of the mass ratio S = mi/m2. The critical mass ratios 6c,i and 6c,2 for L = 1, 2 can 
be read off from the positions where the scaling factor diverges (indicated by the vertical 
dotted lines). For L — )■ oo, the critical mass ratio approaches zero. For 5 — ?■ 0, all scaling 
factors approach unity corresponding to the limit sl — ?■ oo. 

Another important observable is the ratio a*/|a_| that compares the values of the scat- 
tering length a^, and a_ at which Efimov trimers cross the atom-dimer and three-particle 
thresholds, respectively. This ratio can be measured experimentally if at least one resonance 
feature is seen for negative scattering length and one in the atom-dimer system. The calcu- 
lated ratios for L = 1 are shown in Fig. [2] for following one Efimov state and for comparing 
neighbouring states as a function of the mass ratio 6 = mi/m2. 



IV. ANALYTICAL RESULTS FOR ATOM-DIMER SCATTERING 



The scattering of atoms and dimers can be directly related to the STM equation 
equal incoming and outgoing momenta. The scattering amplitude is given by 

k'^^ Ua 

fdk) = , or^i „„. , 7^^ JT^TTT = ^^^k, k; E) , 



^^2L+l 



27r 



with 



(13) 



k^^+^COi5AD,L{k) 

where the energy E and center-of-mass momentum k are related hy E = k^/(2fiA£,)—Bd. The 
dimers are built of two unlike particles and their binding energy is given by Bd = l/(2/xa^) . 
The dimer breakup threshold at E = corresponds to k = \/ jj-ad/ 1^ /ff = ^d- It is useful 
to define an energy dependent scattering length, 
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FIG. 2. The ratio a*/|a_| for following one Efiinov state and for comparing neighbouring states as 
a function of the mass ratio 5 = mi/m2 in the case L = 1. 



Note that this quantity does not have the dimension of a length for L > 0. We can now 
calculate the elastic atom-dimer scattering cross section in the Lth partial wave 



afUk) = {2L + l)A^\h{k)\\ 

The total cross section can be obtained with the help of the optical theorem 

47r, 



a 



(tot) 
AD,L 



(k) = (2L 



D—lmhik) 



(15) 
(16) 



and the inelastic cross section ctad ii^) ^y subtracting Eqs. flT6l) and flTSj) . 

At -E = 0, it is also possible to deduce these quantities analytically using the methods 
of Section IV of Ref. 16|. If the Efimov effect is present, the S-matrix element for elastic 
atom-dimer scattering can be written as 

-e^''"^ cosh(7rsL + isi log(a/ao*) - 7]*)/ cosh(7rsL - isi log(a/ao*) + V*) > (17) 



St 



where ai is a real number and oq* determines the position of the minima in the three-body 
recombination rate for positive scattering length and in the elastic atom-dimer cross section. 
Hence, the elastic cross section can be expressed as 



,(el) 



aTi,{E = 0) = {2L 



^)^\Sl 



''D 



(2L + 1) Ana' 



6{6 + 2) sinh^(7rsL) {sinh^(?7j + sin^[sL log(a/ao*)]} 



{6 + ly sinh^(7rsL + t]^) + cos^[sl log(a/ao*)] 
where kna = {S + 1)/ ^^6(6 + 2) was used. The inelastic cross section is given by 



(18) 



a 



(inel) 






{2L + l)7ia' 



S{6 + 2) 



sinh(27rsL) sinh(2r7#) 



{S + 1)2 sinh^(7rsL + 7]^) + cos'^[sl log(a/ao*)] 



(19) 



Atom-dimer relaxation is the process where an atom and a shallow dimer collide and an 

energetic deep dimer and atom are ejected. It is one of the main loss processes in mixtures 

of atoms and dimers. The atom-dimer relaxation rate constant /3 is defined by the rate 

equation 

d d . . 

-jl'^D = -ji'^A = -pnoUA , (20) 

dt dt 

where ua/d denotes the number densities of atoms and dimers, respectively. The relaxation 
rate is directly related to the inelastic scattering cross section, 

(3l{E) = ^a^-^liE) . (21) 

f^AD 

At E = 0, we therefore obtain 



R (j^-n^-(OT ^^wy-M±3. sinh(27rg^)sinh(2r7,) ha 

/i4i. - U) - (2L + IJTT ^^^^^, sinh\ns, + r^^)+cos^[s,log{a/a,^)]m,- ^''^ 

We note that our analytical results for the inelastic atom-dimer rate at zero energy are in 
agreement with the general scaling law derived by D'Incao and Esry |28|. At the atom-dimer 
threshold, E = —B^, only the S-wave contribution survives. This case was studied in detail 



in Ref. [16[. For all other angular momenta L > 0, the atom-dimer scattering amplitude, 
Eq. (TT3|l . vanishes for fc — ;■ 0. 

Another important process in cold atom experiments is three-body recombination. This 
process can happen in a mixture if not all atoms of species 1 are bound in dimers. The atom 
of species 1 and an atom of species 2 form a dimer, shallow or deep, and another atom of 
species 2 balances energy and momentum. Typically, all three atoms are lost if this process 
occurs in a trap. For angular momenta L > 0, however, the three-body recombination rate 
vanishes at -E = 0. 

In the following, we present numerical results for the atom-dimer observables discussed 
above focussing on the experimentally most relevant case L = 1. 



V. NUMERICAL RESULTS 

A. Atom-dimer observables without Efimov effect 

In this subsection, we only consider the case ?7* = which corresponds to no deeply-bound 
dimers. As a consequence, the inelastic cross section vanishes below the dimer breakup 
threshold. In the P-wave channel, there is no Efimov effect for bosons. The total (=elastic) 
atom-dimer scattering cross section for ii^ = (at dimer breakup) is shown as the solid line 
in Fig. [31 Note that for increasing 6~^, the cross section does not tend monotonically to zero 
but rather oscillates with diminishing amplitude. This oscillation is not due to the crossing 
of bound three-body states with the atom-dimer threshold as three-body bound states are 
not present in this system (see below for the fermionic case). Using Eq. (TTSjl . the elastic 
cross section at threshold can be written as 

^fU^n) = (2L + 1)4™^ ^^ilTi]l si^' ^AoAko) • (23) 
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FIG. 3. The elastic atoin-diiner scattering cross sections for bosons (b) and fermions (f) in the 
P-wave channel and for fermions in the S-wave channel at £^ = as a function of the mass ratio 
6~^. The vertical dotted line indicates the critical mass ratio for fermions in the P-wave channel, 
S'l = 13.61. 



The observed oscillation then implies a monotonic dependence of the elastic phase shift at 



threshold on 5 ^ (29|] with SAD,L{kD) = for 5 ^ = 0. For decreasing energy, the amplitude 
of the cross section gets larger and the peaks move to larger mass ratios 6~^. 

A similar behavior is observed for fermions in the S-wave channel. The atom-dimer scat- 
tering length can be determined according to the formula aAD,o = — ^n(0, 0; —Bd)fiAD/{'^'^) ■ 



This reproduces the results for the mass dependence found by Petrov [2l| that were confirmed 



in Refs. [30|, l3l|- The corresponding total cross section is shown in Fig. [3] as dash-dotted 
line. Again, the oscillation is not due to bound states and Eq. ( 123|) implies a monotonic 
dependence of the threshold phase shift on 5~^. However, in this case SAD,L{k£)) approaches 
a value slightly below 7r/2 for 6^^ = 0. As in the case of bosons with L = 1, we find that 
the amplitude of the cross section gets larger for decreasing energy and the peaks move to 
larger mass ratios 6~^. 

For fermions in the P-wave channel, the Efimov effect only comes into play for mass ratios 
5~^ > 13.61. For the region without the Efimov effect, many observables have already been 



32l-l36j. Kartavtsev and Malykh found one three-body bound state for the 



calculated [2 

range 8.17260 < S'^ < 12.91743 and two three-body bound states for 12.91743 < S'^ < 
13.6069657 (29|. They call these states universal, as their binding energies only depend on 
the dimer binding energy, or equivalently, the^scattering length. The occurrence of these 

They also demonstrated the divergence of 



states was recently confirmed by Endo et al. [36 



the atom-dimer scattering length at the mass ratios where the universal trimer states appear 
and how similar behaviour occurs for higher angular momenta. We have confirmed these 
results. Kartavtsev and Malykh also found that close to the critical mass ratio, the energies 

of the universal states follow a square-root dependence |29|], E — E, 



oc 



^cTi - 5' 



An 



investigation of the behavior of Efimov states for 5 ^ slightly above 5^ \ would be interesting. 



The behavior of the energies must be non-analytic in 5~^ — 5~l, since Efimov states can be 
shifted to any desired energy by adjusting the three-body parameter. However, such a study 
is beyond the scope of the present investigation as our numerical calculations converge only 
slowly close to the critical mass ratio. 

In Ref. [29|, atom-dimer elastic scattering at the dimer breakup threshold was also cal- 
culated. We show this process as the dashed line in Fig. [3] in comparison to the results for 
bosons with L = 1 and for fermions with L = 0. In order to elucidate the physics of the 




5 =m2/m^ 
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FIG. 4. The elastic cross section a^^ ^ in units of a^ versus the mass ratio 6~^ at the energies 
E = 0, —0.91Bii, and — 0.99751?^ (solid, dashed, dash-dotted lines). The dash-double-dotted lines 
indicate the position of the two peaks for varying energy. 



two peaks in the elastic cross section, we show our results for a^^ ^ in Fig. H] as a function 
of 6"^ for three energies, E = 0,-0.915^, and — 0.99755^, i.e., at and below the dimer 
breakup threshold. The dash-double-dotted lines show how the peak positions move with 
varying energy. While the first peak moves monotonically to larger values of S~^ as the en- 
ergy is decreased, the second peak shows a more complicated behavior. Initially, it moves to 
smaller values and reaches a minimum S'^l^^ ~ 12.3 for E / Bd ~ —0.938, before it moves back 
to larger values of 5~^ as the atom-dimer threshold is approached. Our results demonstrate 
that the two-peak structure is indeed due to the presence of the two universal three-body 
bound states. The positions of the two peaks move from 5"^ = 5.63 and 13.31 at £" = to 
two sharp, function like peaks at 6~^ = 8.17 and 12.9 for E = —0.99995^, which are the 
critical values for the occurrence of the three-body bound states. 



B. Atom-dimer observables with Efimov eff'ect 

In the presence of the Efimov effect, the observables do not only depend on the mass ratio 
and the energy but they also depend log-periodically on the scattering length. From now on, 
we focus on the case of fermions with L = 1. We omit the additional subscript 1 indicating 
the P-wave channel for notational simplicity. We find the energy- dependent atom-dimer 



10 



scattering length for 6 ^ > 13.61 calculated with Eq. (fT^ to be very well approximated by 
the formula 



aAD{5,k,a) 



1 



k^ cotdADi^i k, a) 



ci{6, k) + C2((5, k) cot \si log(a/a*) + i?7*] \a^ . (24) 



We show ci/2(5, k) in Fig. [5] for the mass ratios 5 = 0.03, 0.04, and 0.06 as functions of the 
momentum in units of the breakup momentum ko from 0.001 k^ up to the dimer breakup 
threshold. Interestingly, for all considered mass ratios, Ci{6,k) = —8.7 ±0.2 for fc — >■ 0. The 
coefficient €2(6, k) also approaches constant values in this limit but this limit depends on the 
mass ratio. The smaller 6, the larger this approached value. 
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FIG. 5. Coefficients ci{6,k) (dark/blue) and 02(6, k) (light/green) as a function of k/kfj for mass 
ratios 5 = 0.03, 0.04, and 0.06 indicated by the solid, dashed, and dash-dotted lines, respectively. 



At the dimer breakup threshold, we find 



C2{S,k 



D 



k~^ 



and 



ci{6, k 



D 







(25) 



to be very good approximations for 5 < 0.06. This behavior is similar to the case of spinless 
bosons |23|. From Eq. (fT7|) . we can deduce the atom-dimer scattering phase shift Sad at 
£■ = 0. For exp(27rsi) 3> exp(±2?7*), the expression simplifies to 



Sad = 0-1 + 51 ln(a/ao*) + «//* 



(26) 



and the constraints ( l25l) follow from Eq. ( 124|) . Since si approaches as 5 — > 5c,i = 0.07349 
from below, this approximation is invalid at larger mass ratios. In this case, we find that 
Ci(5, ki:,) tends to slightly larger and C2((5, ko) to slightly smaller values. 

The elastic and inelastic atom-dimer scattering cross sections show the typical log-periodic 
dependence on the scattering length. For general momenta, they can be approximated with 
Eqs. (TT3|) - (fT6| [2^ and the appropriate coefficients from Fig. [51 At i? = 0, we can compare 
our numerical calculation to the analytical formulae in Eqs. (IT8|) and (fT9|) . For oad v^d) 
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FIG. 6. Elastic (solid line) and inelastic (dashed line) atoni-dimer scattering cross sections in units 
of a^ for 5 = 0.04, i]^ = 0.1, and E = -0.99i3d versus a|Ac|. 



which shows only a weak dependence on a, we find generally good agreement. For the elastic 
cross section a^l,{k]:,), we find very good agreement in the region 5 < 0.06. For larger 5, the 
numerical calculation becomes difficult because of the large scaling factor. 

As an example, we show the elastic and inelastic cross sections for 5 = 0.04, 77* = 0.1, 
and E = —0.99Bd in Fig. El This mass ratio S roughly corresponds to the mixtures ^Li- 
i7i/i73y]3^ The cross sections show the typical log-periodic dependence on the scattering 
length. The values of the maximal and minimal cross section depend strongly on the energy 
and vary over several orders of magnitude. To demonstrate this dependence, we show the 
maximal and minimal values of a^£, as a function of the center-of-mass momentum k/kjj 
for 5 = 0.03, 0.06 and r/* = 0.2 in Fig. E 

The atom-dimer relaxation rate /3 which can be measured in cold atom experiments is 
determined by the inelastic cross section via Eq. fl2Tl) . In the case of P- waves, /3 vanishes at 
the atom-dimer threshold E = —Bd- In Fig. [HI we show /3 for 5 = 0.04 and 0.06, rj^, = 0.1, 
above threshold for E = -0.99995^, -0.99845^, -0.965^, and 0. We also give the positions 
of the peak in /3 for varying energy. The relaxation rate /3 shows a strong energy dependence 
as well. Starting from being zero at the atom-dimer threshold, it develops resonant log- 
periodic structures for larger energies which become less pronounced as the dimer breakup 
threshold is approached. The peak position is also strongly energy-dependent and varies by 
a factor of four for 6 = 0.04 and by a factor of three for 6 = 0.06. For other mass ratios the 
qualitative behavior is very similar. 



VI. SUMMARY AND OUTLOOK 

In this paper, we have investigated the Efimov effect for heteronuclear systems of two 
identical particles and a third distinguishable particle in higher partial waves. The unlike 
particles were assumed to have resonant S-wave interactions while the interaction between 
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FIG. 8. Dimer relaxation rate /? for 6 = 0.04 and 0.06 in units of ha/mi for r]^ = 0.1 versus 
a|Ac|. The double-dash-dotted, dash-double-dotted, dash-dotted, and solid lines correspond to the 
energies E = —0.99995^, —0.99841?^, — 0.96i3d, and 0, respectively. The dashed lines show how 
the peak position shifts with the energy. 



like particles was neglected. For even (odd) angular momentum L, the two identical particles 
must be bosons (fermions) for the Efimov effect to occur 15|]. Using an effective field theory 
framework, we have derived a generalized STM equation which describes the off-shell atom- 
dimer scattering amplitude in the total angular momentum channel L. All three-body 
observables can be extracted from this amplitude taken in appropriate kinematics. 

We have derived a transcendental equation for the preferred scaling factor exp(7r/si) 
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for arbitrary L as a function of the mass ratio S. The numerical results agree well with a 
previously derived equation using hypergeometric functions derived by Nielsen and cowork- 
ers [19|. For the experimentally most relevant case of the P-wave Efimov effect, we have 
predicted the ratio of the scattering lengths where Efimov states cross the atom-dimer and 
three- atom thresholds, a*/|a„|. This ratio is independent of the three-body parameter and 
can be measured in experiment. For the S-wave case, Barontini et al. [17| have measured 
the value a^,/|a_| = 2.7 in a K-Rb mixture. The universal prediction for this system is 
a*/|a_| = 0.52. The Efimov features in the experiment by Barontini et al. are not deeply 
in the universal region and thus finite range corrections are likely important for this exper- 
iment. The leading finite range correction is given by the effective range correction, but for 
further improvement the van der Waals tail of the interaction has to be included explicitly. 
Moreover, there is a finite energy shift of the atom-dimer rescattering resonance used in 



Ref. [17| to extract a* that needs to be taken into account [16|. For P- waves, no experiment 
has been carried out to date. 

The measurement of atom loss rates has played a key role for the observation of the S-wave 
Efimov effect in cold atoms [l| and the P-wave Efimov effect in a Bose-Fermi mixture could 
be detected in an analogous way. We have derived analytical expressions for the elastic and 
inelastic atom-dimer cross sections as well as the atom-dimer relaxation rate for arbitrary 
angular momentum L at the dimer breakup threshold. For energies below this threshold, 
we have laid out a framework to calculate these quantities numerically. 

Using this framework, we have explicitly calculated the atom-dimer scattering cross sec- 
tions for Bd < E < in low angular momentum channels without the Efimov effect, i.e., 
bosons in the P-wave and fermions in the S-Wave channel. Furthermore, we have calculated 
the cross section for fermions in the P-wave channel below the critical mass ratio 6~l. The 
cross section shows two peaks due to the appearance of two non- Efimov three-body bound 



states |29|. We have calculated the position of these peaks as a function of the collision 
energy E. 

Focussing on the P-wave fermionic channel above the critical mass ratio, we have nu- 
merically calculated the atom-dimer cross section up to the dimer breakup threshold. The 
cross sections show the typical log-periodic dependence on the scattering length. The max- 
imal and minimal cross section values depend strongly on the energy and vary over several 
orders of magnitude. At the atom-dimer threshold, we found good agreement with our 
analytical results. The atom-dimer cross section below the dimer breakup threshold can 
be parametrized by two universal functions Ci{6,k) and 02(6, k). We have calculated these 
functions for several mass ratios and derived simple analytical expressions for their values at 
the dimer breakup threshold if exp(27rsi) ^ exp(±2?7^,) is satisfied. Finally, we have numeri- 
cally calculated the atom-dimer relaxation rate /3 as a function of the three-body parameter, 
mass ratio and energy. As for the bosonic case [16| , the position of the relaxation maxima 
is strongly energy dependent and not a monotonic function of energy. 

In summary, our calculation provides a basis for interpreting experimental results on 
the Efimov effect in higher partial waves. Due to the fermionic nature of the dimers, the 
preparation of the required atom-dimer mixture for the P-wave case should be feasible. A 
few experimental groups already study heteronuclear mixtures of interest to this work, e.g. 



various Yb-Li mixtures in the groups of Takahashi [37[ and Gupta [38[. Other groups inves- 
tigate heavy species that could be mixed with Li, e.g. fermionic Sr in Grimm's group 39 . 
A natural extension of our work would be to calculate the three-body recombination rate 
for energies away from the dimer breakup threshold. This could in principle be done using 
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the methods of Refs. mi 41 
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